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ABSTRACT 

The fraction of ionized gas in the Warm Hot Intergalactic Medium induces tempera- 
ture anisotropies on the Cosmic Microwave Background similar to those of clusters of 
galaxies. The Sunyaev-Zel'dovich anisotropies due to these low density, weakly non- 
linear, baryon filaments can not be distinguished from that of clusters using frequency 
information, but they can be separated since their angular scales are very different. To 
determine the relative contribution of the WHIM SZ signal to the radiation power spec- 
trum of temperature anisotropies, we explore the parameter space of the concordance 
A-Cold Dark Matter model using Monte Carlo Markov Chains and the Wilkinson 
Microwave Anisotropy Probe 7yr and South Pole Telescope data. We find marginal 
evidence of a contribution by diffuse gas, with amplitudes of Awhim = 10 — 20 /iK 2 , 
but the results are also compatible with a null contribution from the WHIM, allow- 
ing to set an upper limit of Awhim < 43 /j,K 2 (95.4% C.L.). The signal produced 
by galaxy clusters remains at Acl = 4.5 /iK 2 , a value similar to what is obtained 
when no WHIM is included. From the measured WHIM amplitude we constrain the 
temperature-density phase diagram of the diffuse gas, and find it to be compatible 
with numerical simulations. The corresponding baryon fraction in the WHIM varies 
from 0.43 to 0.47, depending on model parameters. The forthcoming Planck data could 
set tighter constraints on the temperature-density relation. 

Key words: Cosmology: cosmic microwave background - cosmological parameters; 
Cosmology - theory; Cosmology - observations 



1 INTRODUCTION to detect the WHIM has concentrated in the identification 

of X-ray absorbers, and around half of the WHIM baryons 
At low redsmfts, close to hall or the baryons in the Universe ,,, , ., ,.,-,,, , , ,. ,, , , . ,, 

could have been identified through this method, leaving the 
have yet to be identified (Fukugita & Peebles, 2004 . Numer- . . . . . . ,. , „_ . ., n(w ,„. ., . 

. . . ' remaining missing baryon fraction at ~ 29 ± 13% (ahull ct 



ical simulations suggest the existence of a Warm-Hot Inter 
galactic Medium (WHIM) phase of mildly non-linear struc 



al. 2012). However, many of these detections at z > are 

either of low statistical significance or controversial. For in- 
tures (Cen & Ostnker, 1999, 2006; Dave et al. 2004; Smith et ,, „ j * *• u -m- * i ^ nn ^ , 

, „„-\ ™ „ „ . , \ ^ » , „„..„n ■. , „ , ^ stance, the X-ray detections by Nicastro et al. 2005 have 

al. 2011; Shull, Smith & Danforth, 2012). Most of the efforts 4 , ' „ JU , . ' ,. , v , \ ,'„„. 

' ' not been conhrmed by later studies (Kaastra et al. 2006, 

Yao et al. 2012). As an observational alternative, Atrio- 

* E-mail- rgs@iac es Barandela & Miicket (2006) suggested that the WHIM could 
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be detected through the temperature anisotropies generated 
on the Cosmic Microwave Background (CMB). Due to ther- 
mal and kinetic motions of ionized gas, Compton scattering 
of CMB photons by free electrons induces secondary tem- 
perature anisotropies that were first described by Sunyaev 
& Zeldovich (1970, 1972, hereafter SZ, Birkinshaw, 1999). 
The thermal SZ (TSZ) temperature anisotropies have a dis- 
tinctive frequency dependence, different from other fore- 
grounds, and have been detected in the direction of many 
known clusters (see e.g. Planck Collaboration, 2011). The 
detection of the TSZ signal associated with the less-dense 
WHIM inter-cluster filaments is rather more challenging ow- 
ing to its small amplitude. Hernandez-Monteagudo, Genova- 
Santos & Atrio-Barandcla (2004) attempted to measure this 
signal though cross-correlations between CMB maps and 
templates of the density field constructed from galaxy cat- 
alogues (see also Suarez- Velasquez et al. 2013a). Genova- 
Santos et al. (2008) detected a temperature decrement to- 
wards the Corona Borealis supercluster, in a position with 
no known clusters but with an overpopulation of galaxies 
(Padilla- Torres et al. 2009), but could not confirm it was 
associated to WHIM gas. More recently, an unambiguous 
measurement of a hot and diffuse gas component outside 
the virial regions of clusters was obtained by combining X- 
ray and CMB data from the ROSAT and Planck satellites 
on the cluster pair A399-A401 (Planck Collaboration, 2012). 
The other SZ component, the KSZ effect, is more difficult 
to detect even towards galaxy clusters since it has the same 
frequency dependence as the intrinsic CMB signal. Only the 
average KSZ effect due to large samples of clusters (Kash- 
linsky et al. 2008, 2010) or galaxies (Hand et al. 2012) has 
been detected. 

Frequency information can not be used to distinguish 
the TSZ anisotropy due to clusters and WHIM filaments, but 
they could be separated using the different scale of their re- 
spective anisotropies. The cluster TSZ contribution is max- 
imum at £ max ~ 3000 (Atrio-Barandela & Miicket, 1999) 
while the thermal and kinematic SZ due to the WHIM is 
in the range £ m ax — 200 — 500, depending on model pa- 
rameters (Atrio-Barandela & Miicket 2006, Atrio-Barandela, 
Miicket & Genova-Santos 2008, Suarez- Velasquez, Miicket & 
Atrio-Barandela 2013b). To successfully separate the WHIM 
from the cluster SZ contribution, the latter needs to be 
well characterized. The analysis of the Wilkinson Microwave 
Anisotropy Probe (WMAP) 7 year data confirmed that the 
SZ power spectrum due to the unresolved cluster popula- 
tion was a factor 0.45 ± 0.06 the theoretical predictions and 
estimates based on numerical simulations (Komatsu et al 
2011). This lower than expected amplitude agreed with the 
independent results from the South Pole Telescope (SPT, 
Keisler et al 2011), suggesting that the current models of 
the intracluster medium overestimates the gas pressure rel- 
ative to X-ray and CMB observations (Komatsu et al. 2011). 
An alternative explanation relies on the fact that the TSZ 
power spectrum depends on the amplitude of the matter 
density fluctuations, as, as Cf h tx (erg) 7 (Komatsu & Ki- 
tayama 1999). If the measured value as = 0.801 ± 0.030 



(Larson et al. 2011) were overestimated, so it would be the 
theoretical prediction for the TSZ amplitude. A lower as 
would also result in a lower amplitude of the primordial 
CMB anisotropies at low £, which then could accommodate 
a WHIM component to compensate for this power deficit. 
In Genova-Santos (2009) we showed that a 2-3% variation of 
the concordance ACDM model value of <T8 allowed a contri- 
bution of the WHIM of amplitude £(£ + 1)Ci/2-k ~ 100 fiK 2 
at £ ~ 200. 

In the present paper we first analyze the combined 
WMAP 7yr and SPT data, the best available data set to 
date, constraining the power spectrum out to £ — 3000, 
in order to search for any WHIM SZ contribution to the 
CMB radiation power spectrum. Due to the difference in 
scale, the WHIM and cluster anisotropies can be separated 
through their power spectrum. Since WMAP probes scales 
below £ < 1000 it is mostly sensitive to the WHIM. SPT 
data, that probes £ ^ 600, can constrain better the cluster 
component. By combining both data sets, we could separate 
both contributions to TSZ anisotropies. Second, we study 
if the WHIM component brings the cluster contribution in 
agreement with the theoretical and numerical expectations. 
In Sec. 2 we describe the WHIM model and its parameters; 
in Sec. 3 we describe the power spectrum analysis imple- 
mented in Monte Carlo Markov Chains (MCMC); in Sec. 4 
we discuss our main results and in Sec. 5 we summarize our 
conclusions. 



2 TEMPERATURE ANISOTROPIES 
GENERATED BY THE WARM HOT 
INTERGALACTIC MEDIUM 

In our model, we define the WHIM as shock-heated inter- 
galactic gas at density contrasts in the range Sb — [1, 100] 
and temperatures 10 5 K< T < 10 7 K. We take a rather low 
upper limit for Sb in order to clearly distinguish the WHIM 
in the less dense inter-cluster filaments that form the cos- 
mic web from the gas surrounding clusters of galaxies. The 
WHIM phase of our model encompasses the fraction of the 
missing baryons with Sb ^ 100. Due to the low upper limit 
of 8b chosen by us, the contributions of the WHIM to the 
CMB anisotropy spectrum in our model must be considered 
as a lower limit to the true anisotropy. 

To compute the TSZ power spectrum of this WHIM 
phase, we assume that at any given point the probability of 
having a filament with a baryon density contrast Sb > 1 is 
given by the log-normal probability density function (PDF, 
Atrio-Barandela & Miicket 2006). The log-normal PDF was 
introduced by Coles & Jones (1991) to describe the non- 
linear distribution of matter in the Universe. It has been 
applied to the column density distribution for neutral hy- 
drogen in the intergalactic medium (Choudhury et al. 2001) 
and it has been found to describe very well the matter statis- 
tics at scales larger than 7fe^ 1 Mpc (Kitaura et al. 2009). In 
the log-normal model, the power spectrum of the baryon 
density inhomogeneities follows the dark matter distribu- 
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Figure 1. Different fits of the phase diagram of the hydrody- 
namical simulations of Kang et al (2005). From bottom to top 
the dotted, dashed, dot-dashed and solid lines correspond to val- 
ues 1,2,3 and 4 of the a parameter given in cq. (1). 



tion above a linear scale Lq. Below this cut-off scale, baryon 
density perturbations are smoothed out due to physical ef- 
fects like Jeans dissipation or shock heating. At the physical 
conditions present in the WHIM, shock heating is the dom- 
inating process and Lo(z) is the scale at which the linear 
peculiar velocity v p (x, z) is equal or larger than the sound 
speed c s (z) of the baryon fluid at redshift z (for details, see 
Suarez- Velasquez et al. 2013b). The sound speed is given by 
the mean temperature of the Intergalactic Medium (IGM) 
as c s — (kBTiGM(z)/m p ) 1 ' 2 , where m p is the proton mass. 
At every redshift, the IGM temperature varies with loca- 
tion, but the average value, Tigm, determines the small- 
est baryon density perturbation that survives shock heat- 
ing. The evolution of the IGM temperature is mostly deter- 
mined by the UV background. The mean IGM temperature 
variation with redshift is small and we approximate it by 
log 10 (Ti GM /10 3 K) = (A+0.1(l+z)) (Theunsetal. 2002). At 
redshifts z ^ 3 the temperature varies in the range Tigm = 
10 3 ' 6 - 10 40 K (Tittley & Meiksin, 2007); consequently, 
A = 0.5 - 0.9. At redshifts z > 3 the WHIM does not gener- 
ate significant temperature anisotropies and its contribution 
can be ignored. When Tigm increases, more baryon fluctu- 
ations are erased and the size of baryon filaments increases, 
subtending a larger angle and giving rise to CMB tempera- 
ture anisotropies at larger angular scales. The corresponding 
scale at z = is L ~ 1.7(Ti G M/10 40 if) 1/2 /i _1 Mpc. 

To compute the TSZ temperature anisotropies, it is nec- 
essary to determine the number of free electrons in each 
filament and their average temperature. The number den- 
sity of electrons n e can be obtained by assuming equilib- 
rium between ionization (collisional and photo-ionization) 
and recombination. At temperatures T = 10 5 — 10 7 K and 
density contrasts 8b ^ 100, the gas can be considered fully 
ionized. Whereas in Atrio-Barandela & Miicket (2006) we 
used a polytropic equation of state, in Suarez- Velasquez et 



al. (2013b) we used fits of the phase diagrams obtained in 
various hydrodynamical simulations to take into account the 
effect of shock heating in the evolution of the WHIM. We 
used 



*»gg) 



log 10 (4 + ^+V«) ' 



(1) 



since this numerical fit reproduces well the phase diagram of 
Kang et al. (2005). In this expression, £ = n e /fiB — 8b + 1 
is the electron density in units of the mean baryon density. 
The phase diagrams derived from the simulations are not 
simple linear relations of the type T — T(£), but have a 
large scatter. The a parameter of the previous equation is 
introduced to model this uncertaintly. We consider an inter- 
val wide enough, from a — 1 to a — 4, in order to properly 
cover all the possible variations found in the simulations. 
In Fig. 1 we represent the resulting equations of state for 
q = (1, 2, 3, 4) in decreasing order from top to bottom. 

In our model, the correlation function of the spatial 
variations of the electron pressure along two lines of sight 
separated by an angle 9 is 



C(0) = 



(S(xi,Zl)S(X2, Z2))dzidZ2- 



(2) 



In this expression S — G{y)(k-B<JT /m e c?)n e T e {dU/dzi) with 
i = (1,2). Integrations are carried out along the line of 
sight I. For each direction, the TSZ WHIM temperature 
anisotropy is AT = y c G{v), with y c — feffT / 'm e c 2 J T e n e dl 
the integrated Comptonization parameter, n e the electron 
density, T e the electron temperature, m e c 2 the electron anni- 
hilation temperature, fee the Boltzmann constant, cjt Thom- 
son cross section, dl the line element, dz its corresponding 
redshift interval and G(v) = (x coth(cc/2) —4) the frequency 
dependence of the TSZ effect with x — hv/kTo the reduced 
frequency and To the CMB temperature. Expanding eq. (2) 
in Legendre polynomials gives the WHIM contribution to 
the power spectrum of CMB temperature anisotropies 



Cii = 2tt 



C(6)Pt(cos6)dcosl 



(3) 



where Pi is the Legendre polynomial of multipole I. 

The radiation power spectrum depends on the underly- 
ing cosmological model and the parameters describing the 
physical state of the WHIM, as shown in Fig. 2. The ACDM 
cosmological parameters have been fixed to the best fit val- 
ues of WMAP 7yr data. In Fig. 2a we represent the spectra 
for different fits of the equation of state. The IGM temper- 
ature, that parametrizes the cut-off length of the baryon 
power spectrum Lq, is Tigm = 10 3 ' 6 K, while the amplitude 
of matter fluctuations is as — 0.7. From bottom to top the 
fit parameter of the equation of state is a — 1,2,3,4. At 
densities £ Js 2 the temperature of the free electrons in the 
WHIM increases with increasing a and the resulting tem- 
perature anisotropies also increase. The filled circles show 
the position of the maxima, that in this case vary in the 
range I = max [270,350]. The variation is smaller for larger 
as. These maxima correspond to angular scales in the range 
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Figure 2. Radiation power spectra for different WHIM model parameters, (a) Variation with the a parameter given in cq. 1. From 
bottom to top the dot-dashed, dashed, dotted and solid lines correspond respectively to a = 1, 2, 3, 4, with Tigm = 10 3 ' 6 K and <X8=0.7. 
All other parameters are fixed to the values of the concordance ACDM model. The maxima, indicated by filled circles, are in the range 
£ = [270, 350]. (b) Variation with the IGM temperature and ag for a = 3. Solid and dashed lines correspond to Tjgm = 10 3 ' 6 , 10 40 K; 
for each pair, the top and bottom lines correspond to <rg = 0.7, 1.0, respectively. 



~ 40 — 50 arcmin, much larger than ~ 7 arcmin, the scale 
of the power spectrum of the unresolved 0.5-2 keV Cosmic 
X-ray Background due to both clusters and diffuse gas (Ron- 
carelli et al. 2012, Cappelluti et al. 2012). The difference in 
scale between the power spectrum of TSZ anisotropies due 
to clusters and due to the WHIM is not in contradiction 
with our result. The X-ray power spectra includes the con- 
tribution of diffuse gas from denser environments (8 ^ 10 3 ); 
excluding those high dense regions shifts the power of the 
WHIM anisotropies scales larger than those of clusters. Fur- 
ther, since the TSZ distortion is proportional to the elec- 
tron pressure, n e T e while X-ray emissivity is proportional 
to nlTl 12 (Birkinshaw 1999), the TSZ anisotropy will be 
more extended than the X-ray emission of the WHIM, as it 
occurs with clusters (Atrio-Barandela et al 2008). 

In Fig. 2b we represent the variation with Tigm and 
erg . The solid and dashed lines correspond to Tigm = 
10 3 ' 6 ,10 40 K, respectively. For each pair, the top and bot- 
tom lines correspond to as = 1.0 and 0.7. The fit parameter 
of the equation of state is fixed at a = 3. This figure shows 
that the TSZ WHIM temperature anisotropies grow with in- 
creasing erg and decreasing Tigm- Lowering the cut-off length 
results in larger anisotropies since more baryon density per- 
turbations at small scales survive shock heating. As a result, 
the average angular size of baryon filaments decreases and 
the power spectrum is shifted to small angular scales. 



3 MARKOV CHAIN MONTE CARLO 
ANALYSIS OF WMAP AND SPT 

The shape of the SZ power spectrum, £(£+ 1)Ci/2tv, is very 
similar for both clusters and the WHIM; it presents a single 
maximum at < m „, determined by the average angular size 
of clusters and baryon filaments. The main difference is the 



location of their maxima. As indicated in the introduction, 
it is this difference on angular scale that can be used to sep- 
arate the contribution of the WHIM from that of clusters. In 
both cases, the overall shape of the spectrum is very weakly 
dependent on other model parameters. For instance, due to 
the uncertainties on the number density and pressure profile 
and on their scaling with cluster mass and redshift, Keisler 
et al. (2011) used a fixed spectral shape, derived from nu- 
merical simulations that contains both the TSZ and KSZ 
components, and fit its amplitude to the data. If the radi- 
ation power spectrum is written as Cf z = AszG{v) 2 f(l), 
then f(£) is kept constant and only Asz is derived from the 
data. To search for any WHIM SZ component in addition to 
the SZ due to clusters we will explore the parameter space 
of the standard ACDM model by fitting the theoretically 
predicted radiation power spectra to WMAP and SPT data 
including WHIM and cluster contributions. We use the COS- 
MOMC package (Lewis & Bridle 2002), which implements a 
MCMC method that performs parameter estimation using 
a Bayesian approach. We modified the Keisler et al. (2011) 1 
version by adding at each step of the chain the precomputed 
WHIM SZ component to the theoretical power spectrum 
computed with CAMB (Lewis et al. 2000). 

We constructed a three-dimensional grid of WHIM 
models for values of as = 0.7 and 1.0, a = 1,2,3,4 and 
log(TiGM/10 3 iO = A + 0.1 with A = [0.3, 1.1] varying 
in units of 0.1. Firstly, the power spectra for any given 
value of these parameters were interpolated between the 
pre-computed spectra corresponding to the nearest adja- 
cent values. The interpolation was linear in A and logarith- 
mic in a, as. For the latter, the spectra scales as (7] VHIM oc 



1 Likelihood code downloaded from 
http://pole.uchicago.edU/public/data/keislerll/#Likelihood 
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<j|' 6 (r2B^) 2 (Suarez- Velasquez et al. 2013b). Even if the in- 
terpolation procedure was accurate (errors below 1%), the 
chains were slowly convergent and degenerate with respect 
to Tigm- For this reason, we removed the interpolation on 
the temperature and fixed Tigm • Fig. 2 illustrates that vary- 
ing a (Fig. 2a) and as (Fig. 2b) changes the amplitude 
but does not significantly modify the shape of the radia- 
tion power spectrum. To see the effect of the location of 
the maximum in our final results, we considered two fixed 
WHIM power spectra with maxima located at £ ma x = 330 
and 210. Neglecting the small variations due to the different 
values of a, this maxima correspond to Tigm = 10 ' and 
10 40 K, respectively. Therefore, similarly as it is done for 
the cluster contribution, in our MCMC's we input these two 
spectra and fit their amplitude (Awhim) to the data. 

We fit the data, through the combined WMAP and 
SPT likelihoods, to a concordance ACDM model, de- 
fined by a spatially flat Universe with cold dark matter 
(CDM), baryons, and a cosmological constant A. Following 
Keisler et al. (2011) we added the combined galaxy clusters 
TSZ+KSZ amplitude (Acl) and two amplitudes associated 
with the contributions from Poisson (Apg lsson ) and clustered 
(Ap S ustcred ) distributed point sources in addition to the am- 
plitude of the WHIM TSZ signal (Awhim). The WHIM and 
cluster power spectra were normalized to unity at £ = 300 
and £ = 3000, respectively. Then, to take into account the 
frequency dependence of the TSZ effect, their amplitudes 
Awhim and Acl were scaled to the average frequency of 
WMAP V and W bands ((G(v) 2 ) = 2.861) and of the three 
SPT bands {{G{v) 2 ) = 1.107). Note that in our analysis 
we are neglecting the less important KSZ component of the 
WHIM. Appart from being fainter, this is justfied because 
the TSZ and KSZ WHIM power spectra have very similar 
shapes, rendering a joint fit of both components very degen- 
erate. 

We computed MCMCs in three different cases: to re- 
produce the analysis of Keisler et al. (2011) one run con- 
tained only the clusters SZ contribution; two other runs 
included both cluster and TSZ WHIM contributions, each 
with a fixed value of the IGM temperature, Tigm = 10 3 ' 6 
and 10 40 K. For each of these three cases, we ran eight inde- 
pendent chains with a total number of ~ 350, 000 samples. 
The Gelman & Rubin (1992) criterion showed that all our 
chains had converged; in all cases the R statistic was well 
below 1.2. For instance, R ss 1.01 for Awhim and « 1.003 
for the other parameters, when the WHIM component was 
included, and R ~ 1.001 when it was not. Keisler et al. 
(2011) used a Gaussian prior in the SZ cluster amplitude, 
A CL = 5.5 ±3.0 ^K 2 , derived by Shirokoff et al. (2011) from 
an analysis of an earlier SPT data release. To determine the 
bias introduced by this prior on the WHIM amplitude, we 
ran our MCMCs with and without prior. No significant dif- 
ferences were found (except, of course, on the best-fit value 
of Acl itself) and henceforth we will quote results with no 
prior. 



4 RESULTS AND DISCUSSION 

In Fig. 3 we represent the marginalized likelihood functions 
of the amplitude of the TSZ anisotropy generated by WHIM 
(Fig. 3a), TSZ and KSZ anisotropy generated by clusters 
(Fig. 3b), physical baryon density (Fig. 3c) and of the ampli- 
tude of matter density perturbations at 8ft _1 Mpc (Fig. 3d). 
For an easier comparison, all likelihoods were normalized 
to unity. As mentioned above, the results presented do not 
include a prior on Acl- Blue and red lines correspond to 
Tigm = 10 3 ' 6 and 10 40 K, while the cyan lines correspond 
the results when no WHIM is included. Adding a WHIM 
component marginally decreases Acl and as- In Fig. 4 we 
plot the two-dimensional likelihoods for different pairs of pa- 
rameters. The figure shows that Awhim is degenerate with 
respect to other parameters, especially with respect to Acl 
and as, and that it remains compatible with zero at the 2- 
a level. The 1-a contours of the Awhim — as plot appear 
slightly curved downwards for increasing Awhim, reflecting 
the (marginal) decrement of as when we include the WHIM 
TSZ anisotropy. Finally, the right panel shows how fie ft 2 
decreases when there is no WHIM contribution while Acl 
remains almost unaffected. The dashed lines show the con- 
tours when no WHIM anisotropy is included. 

In Fig. 5 we compare the measured WMAP 7yr (red 
filled circles) and SPT (blue filled circles) with our best- 
fit models in two cases (a) with a WHIM anisotropy with 
Tigm = 10 3 ' 6 K and (b) without a WHIM component, repre- 
sented by (indiscernible) solid lines. We also plot the differ- 
ent components: primordial CMB (dashed lines), TSZ from 
WHIM (dotted line), SZ from galaxy clusters (dashed-dotted 
lines) and emission from point sources (both Poisson and 
clustered terms; dashed-triple-dotted lines). For visual pur- 
poses, we have scaled the WHIM component to the WMAP 
frequencies and the cluster component to the SPT frequen- 
cies. 

Table 1 shows the parameters that define our best- 
fit cosmological model for our three MCMC (without Acl 
prior). These parameters are derived from the median of 
the marginalized posterior probability density function, and 
their confidence intervals encompass the 68% of the proba- 
bility around those points. The parameters derived when no 
WHIM component is included differ by less than 2% from 
those of Keisler et al. (2011). In the three cases considered 
by us, the amplitude of the SZ signal from clusters or from 
the WHIM is compatible with zero at the 2a level. From the 
marginalized likelihood distributions we derive upper limits, 
at the 95.4% C.L., of Awhim < 43.6 ^iK 2 and < 42.7 ^K 2 , 
for Tigm = 10 36 K and 10 40 K, respectively. Taking the 
most-likeliy Awhim at face value, from the measured values 
of Q.bIi 2 , Hq and as given in the table, we can determine the 
value of a that reproduces the amplitude of Awhim for each 
value of Tigm- By interpolating in our grid of models we find 
a = 2.46, 3.40 for Tigm = 10 3 ' 6 , 10 40 K, respectively. 

For each of these values of a we can compute the baryon 
content residing in the WHIM filaments with overdensities 
in the range [^1,^2]- To this aim, we use the formalism de- 
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Figure 3. Marginalized likelihoods for the WHIM TSZ amplitude at I = 300 (Awhim), for the combined TSZ+KSZ amplitude at 
£ = 3000 from galaxy clusters (Acl)i for the physical baryon density (Cl-gh 2 ), and for the amplitude of matter density fluctuations in 
scales of 8 h _1 Mpc (ag). We plot our results for the two IGM temperatures, Tigm = 10 3,6 K (blue) and 10 40 K (red), as well as when 
no WHIM signal is introduced (cyan). 
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Figure 4. Mean two-dimensional likelihoods for three combination of parameters, for Tigm = 10 3 ' 6 K. Solid lines represent the la and 
2a confidence regions. Dashed lines in the right panel correspond to the case with no WHIM component. 



scribed in Suarez- Velasquez et al. (2013b). The mass frac- 
tion is given by M(£i,£ 2 ) = £ _1 f^ £-F(0> where £ is the 
mean of the distribution and F (£) is the log-normal proba- 
bility distribution function. The integration range must be 
limited to those overdensities where the non-linear evolu- 
tion is well described by our model. As a criteria, we took 
£2 ^ 100 and £1 = £ mo dian + o", with a the standard deviation 
of the log-normal distribution. To include only scales that 



undergo shock heating we also require as necessary condition 
that fi > 2 at all redshifts. Since £1 and Tigm are redshift 
dependent, the fraction of baryons residing in the WHIM fil- 
aments is a decreasing function of z. For a — 2.46, 3.40 and 
Tigm = 10 3 ' 6 , 10 40 K the derived baryon fractions at z = 
are respectively fiB.WHiM/ffe ~ 0.43 and 0.47. Notice that 
even if the value of a is rather uncertain, the baryon frac- 
tion varies little within the allowed parameter space. Our 
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Figure 5. Experimental power spectra measured by WMAP7 (red dots) and SPT (blue dots) compared with our fitted models. Dashed, 
dotted, dashed-dotted and dashed-triple-dotted lines correspond respectively to the primordial CMB, WHIM tSZ, combined tSZ and kSZ 
signals from clusters, and point sources (combing the terms coming from Poisson and clustered distributed sources). Solid lines represent 
the sum of all these components. We compare our results including the WHIM TSZ component (magenta curves) with the case when no 
WHIM component is included (cyan curves). This plot corresponds to the case of Tigm = 10 3 ' 6 K, and shows no significant differences 
with respect to the case of Tigm = 10 40 K. 



result is slightly above the 0.29 ± 0.13 fraction of baryons 
that remain unidentified according to Shull et al. (2012). 
The baryons detected through TSZ anisotropics will include 
baryons on low dense regions identified by other techniques, 
like some of the OVI absorption systems; after correcting the 
double accounting, the fraction of identified baryons would 
be close to unity. 

In our chains we have fixed Tigm; then, including the 
WHIM contribution adds only one parameter, Awhim (that 
is related with a), to the chain. The improvement in the 
quality of the fit is minimal: in the best case A^ 2 — —0.4. 
Therefore the data fails to provide any significant evidence 
of the existence of the WHIM. Not unexpectedly, when in- 
cluding a WHIM anisotropy, Acl decreases by 10%, en- 
hancing the disagreement between the measured cluster SZ 
anisotropics and the theoretical expectation. However, the 
discrepancy is alleviated since as also decreases. In the 
model with Tigm = 10 3 ' 6 K, the lower value of as reduces the 
theoretical SZ signal by 7% compensating the lower value of 
Acl found. To make the cluster SZ signal compatible with 
the theoretical predictions we would need as ~ 0.73, much 
smaller than the measured value. 

The results of Table 1 can be understood taking into ac- 
count the difference in the angular scales probed by WMAP 
and SPT. The information on the radiation power spectrum 
provided by these instruments is complementary with each 
other. In WMAP, at low ^'s errors are dominated by sam- 
pling variance and at high £'s by noise. By contrast, the 
radiation spectrum has been best measured by the SPT 



in the range £ ~ 600 - 3000. Then, while WMAP data is 
mostly sensitive to the range where the WHIM contribu- 
tion is largest, around the first and second acoustic peaks, 
SPT data constrain better the cluster TSZ contribution. 
The first acoustic peak occurs at £ ~ 200, its amplitude 
is 1(1 + l)Ce/2ir ~ 6000 fj,K 2 and scales approximately as 
af. In ACDM if all parameters are held fixed but as de- 
creases by 0.5%, the amplitude of the first acoustic peak 
decreases by ~1%. i.e., there could be a WHIM as high 
as 60 /xK 2 at the maximum. Since WMAP operates at 
Rayleigh- Jeans frequencies, this corresponds to an ampli- 
tude Awhim = 15 — 20 (J.K 2 , compatible with the results 
of Fig. 3. Even though the WHIM contribution falls for 
£ ^ 1000, the tail of the distribution overlaps with clus- 
ter anisotropics and the SPT data would suppress the latter 
(see Fig. 5). 

Since our results are compatible with zero WHIM con- 
tribution, one could interpret them as an indication (1) that 
the data has not enough statistical power to identify the 
WHIM contribution or (2) that the WHIM is not a strong 
contributor to the baryon budget. While the estimated val- 
ues of the total baryon mass fraction in groups and clusters 
are still lower than the latest CMB measurement of the same 
quantity (Giodini et al. 2009), it seems implausible that a 
better modelling of the gas physics could eventually solve the 
WHIM problem. If clusters and groups store more baryons 
than presently believed, the discrepancy between the SZ am- 
plitude of clusters measured by SPT and WMAP with the 
numerical expectations summarized in Komatsu et al (2011) 
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Parameter 


No WHIM 


WHIM T IGM = 10 3 - 6 K 


WHIM T IGM = 10 40 K 


n B h 2 


0.0223 ± 0.0004 


0.0228 ± 0.0005 


0.0226 ± 0.0005 


^CDM^ 2 


0.110 ±0.005 


0.111 ±0.005 


0.111 ±0.005 


1000 


1.041 ±0.002 


1.042 ±0.002 


1.042 ±0.002 


r 


U - US5 ' -0.0070 


0-0864+™ 


0.0865^™ 


n s 


0.964 ±0.011 


0.962 ±0.011 


0.967 ±0.011 


ln(10 10 A s ) 


3.19 ±0.04 


3.18 ±0.04 


3.17 ±0.04 


Y He 


0.2478 ± 0.0002 


0.2480 ± 0.0002 


0.2479 ± 0.0002 


-4\VHIM 




19 l 3fi + 1 3' 37 

±».ou_ 13 34 


19 16+ 12 ' 65 
iy.io_ 12 91 


A C h 


^•15-3.37 


4.58±l;ii 


4.60t|;i 2 


zlPoisson 


20.40iJ|3 


20.23 ± 2.89 


20.32t2.| 7 


A Clustered 
^PS 


t- -,,+2.22 
- ±J --2.28 


c- qo+2.21 
o.uo_ 2 25 


5 m +2/20 

o.uo_2 26 


^8 


0.803 ± 0.024 


0.795 ±0.025 


0.800 ± 0.024 


n A 


0.735 ±0.025 


0.738 ± 0.025 


0.736 ±0.025 


f^M 


0.265 ±0.025 


0.262 ±0.025 


0.264 ±0.025 


Zre 


10.40±1;1£ 


10.33 ±1.18 


10.39l;i i? 


Ho 


71-00+llt 


71 5S+ 2 ' 19 
d.OO_2 20 


71.27 ±2.14 


x 2 


3756.04 


3755.84 


3755.66 



Table 1. Concordance ACDM best fit parameters for each MCMC. We show the results for the three cases that we have considered: no 
WHIM, and a WHIM component with Ttqm temperatures of 10 3,6 and 10 40 K, respectively. Parameters are derived from the median 
of the marginalized posterior probability density functions, while the confidence intervals are calculated from the 68% area around those 
points. The last row shows the \ 2 of these models. The amplitudes AwhiMj ^CLi Apg lsson , ^Clustered are gj vcn m U nits f ^K 2 and the 
Hubble constant Hq in kms -1 Mpc~ 1 . 



would be even more acute. What our results indicate is the 
difficulty to accommodate a WHIM component that brings 
the measured cluster TSZ signal in agreement with the the- 
oretical predictions. The reason of this discrepancy might 
not be the measured value of as or Asz, but the lack of 
understanding of some of the aspects of the cluster physics 
that are introduced in the simulations. 



5 CONCLUSIONS. 

We have explored the parameter space of the concordance 
ACDM model using two SZ contributions: one due to the 
combined KSZ and TSZ effects from the unresolved cluster 
population, the other due to the TSZ from WHIM. We have 
fitted this model to the combined WMAP 7yr and SPT data, 
at present the best data-set publicly available, covering the 
multipole range £ = [2, 3000]. We have found that a WHIM 
component with an amplitude of ~ 20 fiK 2 at £ ~ 300 is 
compatible with the data. This new WHIM component re- 
sults in a 10% decrement in the amplitude of the SZ signal 
ascribed to galaxy clusters, while the cosmological parame- 
ters do not change significantly; QbIi 2 is 2% higher, and as 



decreases by 1%. This new fit does not solve the discrepancy 
between the measured amplitude of the cluster power spec- 
trum and the theoretical expectation; this would require a 
considerably lower value of as- However, since as is tightly 
constrained by the data, this could be an indication of im- 
proper modelling of cluster evolution, gas dynamics, and of 
the resulting pressure profiles, both theoretically and numer- 
ically. Our analysis shows that a WHIM contribution is com- 
patible with WMAP and SPT data and, if this contribution 
is well described by our model, the average properties of the 
WHIM can be measured. Since the cosmological parameters 
(08, $!b and h) are determined from the CMB anisotropy 
data, measuring the WHIM contribution allows to constrain 
the other two parameters, namely, the cut-off length of the 
baryons power spectrum, parametrized by Tigm, and the 
phase diagram parameter a that also determine the ampli- 
tude of the WHIM anisotropy. For Awhim ~ 20fiK 2 the 
Kang et al. (2005) phase diagram parameter is constrained 
to be in the interval a ~ [2.5,3.4], consistent with the esti- 
mates derived from simulations. The fraction of baryons that 
would reside in this phase, between 43% and 47%, would be 
large enough to close the baryon census problem summa- 
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rized in Shull et al. (2012). These are the first constraints 
on the physical state and abundance of the WHIM presented 
to date. 

Planck, with its wide coverage of frequencies and an- 
gular scales, could help to unambiguously measure this 
WHIM component, separating it from the cluster contri- 
bution, and to set more stringent constraints on the WHIM 
physical properties. The power spectrum amplitude scales 
as cf mM oc AwhimG(^) 2 . If Awhim = 10-15 ^K 2 at 
Planck frequencies it will change from zero at 217 GHz to 
~ 60 /xK 2 at 40 GHz (decrement) and 353 GHz (increment). 
This represents a 2% variation across the observed range. If 
foreground contributions, that also change with frequency, 
are subtracted down to this level, then the variation of the 
power spectrum will be an indication of the WHIM TSZ 
contribution. 
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